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Abstract 

Trapped surfaces are studied as inner boundary for the Einstein 
vacuum constraint equations. The trapped surface condition can be 
written as a non hnear boundary condition for these equations. Un- 
der appropriate assumptions, we prove existence and uniqueness of 
solutions in the exterior region for this boundary value problem. We 
also discuss the relevance of this result for the study of black holes 
collisions. 

PACS: 04.20.Dw, 04.20.Ex, 04.70.Bw 

1 Introduction 

A black hole defines a boundary in the spacetime: the boundary of the region 
of no escape to infinity. This boundary is called the event horizon. There exist 
others kinds of spacetime boundaries, for example the ones defined by matter 
sources: the interface between matter and vacuum. In this case the boundary 
is introduced in the sources of Einstein equations by choosing an energy 
density v^ith compact support. For black holes, the boundary is produced by 
the vacuum equations themselves, it depends only on fundamental properties 
of gravity; in this sense it is a more fundamental kind of boundary than the 
matter sources ones. 

In the presence of a boundary, is natural to study boundary conditions 
for the equations and solve only for the exterior region. In the case of black 
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holes this is very desirable since, firstly, the observations are made at infinity 
which is causally disconnected with the black hole region. And, secondly, in 
the black hole region there are singularities which are very difficult to handle 
in the numeric simulations. 

In the context of an initial value formulation, the first step in order to 
understand the black hole boundary value problem is the study of the inter- 
section of this boundary with a spacelike three-dimensional Cauchy hyper- 
surface. That is, to study the black hole boundary value problem for the 
constraint equations. This will be the subject of this article. 

There is no known way to characterize the intersection of the event hori- 
zon with a Cauchy surface in terms of a differential condition which involves 
only fields on the Cauchy surface. Such local characterization is provided 
by the concept of trapped surface. A trapped surface is a compact, two- 
dimensional surface for which the expansions of both sets, outgoing and in- 
going, of future directed null geodesies orthogonal to the surface are negative. 
The relevance of trapped surfaces relies in two important results: i) a trapped 
surface is always inside an event horizon, ii) The development of an initial 
data set which contain a trapped surface will be geodesically incomplete (see 
|26] and (SH])- Because of i) we expect that trapped surfaces should be good 
boundaries for the constraint equations. In this article we will see that in 
fact this is true: trapped surfaces provide inner boundary conditions for the 
constraint equations under appropriate assumptions. 

Associated with trapped surfaces there is the concept of an apparent 
horizon. An apparent horizon is essentially defined as the outermost trapped 
surface. In this article only trapped surfaces will be discussed and not ap- 
parent horizons, since with the techniques used here it will allows to decide 
whether a given trapped surface is in fact the outermost. 

The problem of constructing initial data with several disconnected trapped 
surfaces is relevant for the study of black holes collisions (see for example [TT] 
and [ini)- In almost all data currently used in numerical simulations (the 
exceptions are the numerical studies in [36j and trapped surfaces are 
produced by indirect means, for example introducing a non trivial topology. 
The theorem proved in this article will provide a way of solving only for the 
exterior region, with the appropriate boundary conditions, without imposing 
any symmetry and without introducing any non trivial topology. 

From the point of view of getting initial data for black hole collisions 
this has several advantages. Firstly, no computer resources are wasted in 
the interior region. Secondly, the location of the trapped surfaces is known 
a priori. Finally, the control of the boundary conditions presented here will 
allow one to construct more general classes of black hole initial data than the 
ones studied so far (see the discussion in section [21). 
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The plan of the article is the following. In section [21 we present our result, 
given by theorems 12. II and 12.21 We also discuss the physical interpretation of 
these theorems. In section El we collect some results from the elliptic theory 
that will be used in the proofs. In section |31 and section Owe prove theorem 
12.11 and theorem 12.21 respectivelv. 

2 Main Result 

Let Ck be a finite collection of compact sets in M'^. We define the exterior 
region Vt = ]R^\UfcCfc. An initial data set for the Einstein vacuum equations is 
given by the triple (fi, hab, Kab) where hab is a (positive definite) Riemannian 
metric, and Kab a symmetric tensor field on 17, such that they satisfy the 
constraint equations 

b'kab - Dak = 0, (1) 

R + k^ -kabk^^ = (2) 

on f2; where Da is the covariant derivative with respect to hab-, R is the trace 
of the corresponding Ricci tensor, and K = h°'^Kab- 

The data will be called asymptotically flat if there exists some compact 
set C, with UfcCfc C C, such that VL\C can be mapped by a coordinate 
system diffeomorphically onto the complement of a closed ball in M'^ and 
we have in these coordinates 

— 27T7 

/i,, = (l + — )5,,+0(r-2), (3) 

Ki, = 0(f-2), (4) 

as f = (X]j=i(^"')^)^^^ ~^ ^1 where the constant m is the total mass of the 
initial data. 

The boundaries dCk are assumed to be smooth, two dimensional surfaces 
in {Q, h). Let be the unit normal of dCk, with respect to hab, pointing in 
the outward direction oifl. Let be the unit timelike vector field orthogonal 
to the hypersurface fl with respect to the spacetime metric gab {tH^Qab = 
—1 with our signature convention) The outgoing and ingoing null geodesies 
orthogonal to dCk are given by l'^ = t"' — z/" and fc" = + z/" respectively, the 
corresponding expansions are given by 0+ = Va^" and 6_ = Va^", where 
Va is the connexion with respect to gab, see Fig. Ill We can calculate these 
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expansions in terms of quantities intrinsic to the initial data 

Q_ = k + H- z>V^i^,fe, (5) 
Q+ = k~H- iy'^iy'kab. (6) 

where H = Dai'"'. Note that K = Vat"" is the mean curvature of the three- 
dimensional hypersurface f2 with respect to the spacetime metric gab and the 
normal and H is the mean curvature of the two-dimensional surface dCk 
with respect to the metric hab and the normal z/". Our choice of the definition 
of the extrinsic curvature is given by Kab = hlVch which agree with IHEj but 
is the negative of the choice made in [33] . If dC is an sphere of radius f on a 
canonical Minkowski slice t = const., then we have 0+ = 2/f, 0_ = — 2/f, 
H = -2/f, kat = 0, and 5" = -{d/dfy. 




Figure 1: The exterior region Q with two boundary components dCi and 
dC2. 

The boundary dCk will be called a future trapped surface if 9+ < and 
6_ < on dCk and a future marginally trapped surface if 0+ < and 

e_ < 0. 

We want to find solutions of the constraint equations JfJ-Q which are 
asymptotically flat (i.e; they satisfy ©-(III)) and such that all the boundary 
components dCk (i.e; the whole boundary dfl) are trapped surfaces. In order 
to achieve this, we will reduce the constraint equations to an elliptic boundary 
value problem in which a negative (non-positive) ingoing null expansion 6_ 
can, essentially, be freely prescribed at the boundary dQ. Under further 
restrictions on 9_ it will follows that also 9+ will be negative (non-positive), 
and hence the boundary will be future trapped (future marginally trapped). 
We emphasize that only 9_ (and not 9+) will be free data on the boundary. 

The elliptic reduction will be given by the conformal method. Also, in- 
stead of working with the exterior region f2 it will be more convenient to 
work with its related compactification Q. In the following we describe both 
procedures (see IH], [IHl and the references therein for a description of the 
conformal method, and references [S], [2111, ISI], [IBI for the compactification 
procedure). 
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Figure 2: The compactification fl of the exterior region fl showed in Fig. [H 



Let fl a bounded domain in M^. We will assume that Q is connected, 
however in general it will not be simply connected, each hole in f2 will cor- 
respond to an extra boundary component dCk- Let i G f2 be an arbitrary 
point and set = f2 \ {i}. The point i will represent the infinity of fl. Let 
hab be a Riemannian metric defined on fl. We denote by u"- the unit normal 
to dfl with respect to hab, pointing in the outward direction, see Fig. El We 
will assume that the conformal metric satisfies 

KbeC^{Q\{i})nC\fi). (7) 

We make different assumptions concerning the regularity of hab 8Ltfl\{i} and 
at i. The regularity assumption at i is concerned with the fall off behavior 
of the solution at infinity, it is not related with the inner boundaries dCk. 
The physical motivation for this distinction has been discussed in [T^, [El 
and [12]; in particular assumption Q is weaker that the ones made in these 
references. 

Let K"'^ be a symmetric, trace free tensor with respect to hab, such that 

DaK'''' = on fi, (8) 
and let be a positive solution of 

LhiJ = -iKabK'^'iJ'' on fl, (9) 

o 

where Lh = D'^Da — R/8 and R is the Ricci scalar of the metric hab- If we 
are able to find solutions of equations JHI)-®, then the physical fields {h, K) 
defined by Kb = i^^hab and K"'' = ^p'^^K"^ will satisfy equations {H)-© on 
fl. Since we have assumed = we will obtain K = 

We discuss now boundary conditions for equations (jH])-©- There are two 
kind of boundary conditions. The first one is asymptotic fiatness, we want to 
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ensure that the solution will satisfies the fall ofi' ©-(III)- For this we require 
and 

limr-?/' = 1 at (11) 

where are normal coordinates with respect to the metric hah centered at i 
and r the corresponding radius. 

The second boundary condition is the requirement that 6_ should be 
freely prescribed at the boundary dQ.. In order to write this condition we first 
calculate the null expansions 9+ and 9_ in term of the conformal quantities 

e_ = i)-^ {Au^'DatlJ + HiIj- tfj-^Kahiy'^iy'') , (12) 
e+ = ^p-^ i-Au'^DatP -H^- ^-^Kabiy^iy") , (13) 

where we have used that K = K = 0, the normal vectors are related by 
i/" = ip~'^u°- (that is, p"" is a unit vector with respect to the conformal metric 
hab) and H = DaV'^. From equation (fT2|l we deduce the following boundary 
condition for the conformal factor ip 

Nhtfj = ^^e_ + tp-^Kabiy^u^ on on, (14) 

where Nh = Au^Da + H. 

Equation (fT^ is a boundary condition for the conformal factor only. Re- 
mains to prescribe the boundary condition for the momentum constraint (jH)). 
From (fT^ we see that the simpler way of solving our problem will be to chose 
for the momentum constraint a boundary condition such that the function 
KabV°'v^ can be freely prescribed on dVL. Our first result, concerning the mo- 
mentum constraint, essentially says that this is possible. In order to write 
this theorem we need to introduce some additional concepts. 

The conformal Killing operator Ch, acting on vectors fields w", is defined 

as 

2 

{ChW)ab = '2Di^aWa) " -habD^Wc- (15) 

We say that ^'^ is conformal Killing vector field if it satisfies ChC, = on Q. 
Given a conformal Killing vector field we define the conformal Killing data 
at the point i by 

ka = ]: DaD^eii), S'^ = e\,D'i\i), = ^ (^), a = ^Z^„r(0- (16) 

D 6 

Since Vt is connected, the integrability conditions for conformal Killing fields 
(cf. [nni) entail that these ten conformal Killing data at i determine the field 
uniquely on f2. 



6 



Because of the singular behavior of K"-'' at i given by (fTTH) . it will conve- 
nient to write the results in terms of a further rescaled metric h'^^ for which 
the curvature vanished at i. The metric h'^^^ can be explicitly computed as 
follows. Denote by the open ball with center i and radius r = e > 0, where 
e is chosen small enough such that is a convex normal neighborhood of 
i. We define the cut-off function Xe as a non-negative, smooth function such 
that Xe = 1 in -Be/2 and Xe = in \ i?^. Consider the following smooth 
conformal factor 

= e^^^o with fo = \ x^x^ L,k{t), (17) 

where we have used the value at i of the tensor Lab = Rab — \Rhab- By a 
straightforward calculation we get that the Ricci tensor of the metric 

Kb = ^Xb (18) 

vanishes at the point i. 

Equation ^ is solved using the standard York splitting (cf. flU^) adapted 
to our setting. The free data is a trace free, symmetric, tensor Q"''. Since we 
have two kinds of boundary conditions, Q"''' has a natural decomposition 

Q ~ Qsing ~^ Qregi (^9) 

where Q'^Jl^g and Q'^'lg can be roughly characterized as follows (see section |31 
for details). 

The tensor Q^^^g different from zero only in B^ and pick up the singular 
behavior of A'"* at z, i.e; it will blow up like at i. It contains the linear 
and angular momentum of the data. These physical quantities appear as 
constants and J"^ which partially characterize Ql\ng- In addition to P" 
and J°, there are other constants and A involved in Q'^ing'i these ten 
constants are related with the conformal Killing data (fT6|l as we will see. 

The tensor Q^g^ gives the boundary value of K"-'' at dfl and will not con- 
tribute to the linear and angular momentum. For example we can take Q^lg 
smooth in 0; however a more general behavior at i is allowed, in particular 
Q^lg can blow up like at i. 

With these definitions, we can write our existence theorem for the mo- 
mentum constraint. 

Theorem 2.1. Assume hab satisfies {^j and dVL is smooth. Let Q"''^ the 
symmetric, trace free tensor, given by {I^ where Q^'ig (indQI^^g satisfy M^)- 
J^J) and respectively. 
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i) If the metric hab admits no conformal Killing fields on ^, then there 
exists a unique vector field G C°°(fi \ {i}) n H^{Vl) such that the tensor 
K^^ defined by 

K''' = cul\Q'''-{CH,wr'), (20) 

satisfies the equation DaK"^^ = in Q and K'^'^Va = Qrlg^a on dfl. The 
metric h'^i^ and the conformal factor uq are defined by ( TT^] and ( fTTj) . 

a) If the metric hat admits conformal Killing fields on a vector field 
w"' as specified above exists if and only if the constants P", J", A and 
(partly) characterizing the tensor field Q^ing, satisfy the equation 

P^K + rSa + Aa + {P'L,^{^ + Q'')qa= I Qf^gUbdS, (21) 

Jan 

for any conformal Killing field ^" of hat; where the constants k"-, S"", a and 
q'^ are the conformal Killing data at i for given by 

In both cases i) and ii) K"^ is unique, K"^ G C°°{n\{i}) and K"^ = 0{r-^) 
at i. 

Note that in case ii) K"-'' is unique but w"- is not, we can add to a given 
solution any conformal Killing vector. Condition (|^T]l . which arise in the 
presence of conformal symmetries, is the natural extension to the analog 
condition discussed in [H] and ^Hl where there is no inner boundary and in 
[T7| where there exists matter sources in a compact region. The fact that 
not only Killing vectors but also conformal Killing vectors appear in (^111 is 
a particular feature of maximal slices (i.e; K = K = 0), see the discussion in 

m 

The next result is concerning the Hamiltonian constraint (jH]). In order to 
enunciate it, we need two auxiliary solutions. The first one, denoted by tpQ, 
is a solution of the problem with K"'^ = on f2 (i.e; time symmetry) and 
6+ = 9_ = on dQ. That is, ipo satisfies the following linear boundary 
value problem 

Lh^po = in 1^, (22) 
Nh^o = on on, (23) 
limr-j/'o = 1 at i. (24) 

r— >0 

The second particular solution, ipi.is a solution of the following linear bound- 
ary value problem 

Lhiji = -liT^.i^'^Vo"' in ^, (25) 

o 

Nf,^, = ^^KabV-v' on (26) 
limripi = 1 at (27) 

r— >0 
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Theorem 2.2. Assume that the conformal metric hab satisfies (0) and that 
dVt is smooth. Assume also that R > in D:, H > on and that either 
R or H is not identically zero. Then: 

i) There exist a unique, positive, solution ipo of ([^ - (f^ . and ipo ^ 

c^in\{t}). 

a) Let K"-^ be given by theorem \2.1\ Assume, in addition, that 

Qtg^aVh > on on. (28) 
Then there exist a unique, positive, solution ipi of (f25|l - pTj) . and ipi G 

c^in\{^}). 

Hi) Let 0- < and assume it satisfies the following inequality 

\e-\<ij^'Qf,gUM. (29) 

Then there exist a unique, positive, solution ip of equation ^ with boundary 
conditions (fTT|l and Jill), andip G C°°(VL\{i}) . Moreover, ipo < ip < ipi and 

e+ < e_ < 0. (30) 

Let us discuss the assumptions and the conclusions of this theorem. We 
have assumed that the boundary dfl is smooth, this is done for simplicity and 
because there seems not be any physical reason to study trapped surfaces 
with rough boundary. For rough boundaries (for example with corners or 
cusps) the following proofs are not valid. 

The condition R > and Dai'"' > ensures that the conformal metric 
hab can be conformally rescaled to a metric hab for which the boundary dfl 
is an extremal surface (i.e; H = 0) and R > 0. Then, this condition can 
expressed in terms of the corresponding Yamabe class studied in [TH|. A 
similar condition also appears in the case without inner boundary, see and 
also [21], [ini for a discussion of the same condition using a compactification 
as we use here. 

Simple examples of sets {h, fl) which satisfy these conditions are the fol- 
lowings. Take hat = Sab, where 6ab is the flat metric. Let f2 = 5,.,,, where 
is a ball of radius tq. Then H = 2/ro and R = 0. This example 
has only one boundary component. Note that, in this example, we can not 
take a second ball inside Br^ to construct a region like the one showed in 
Fig. 121 because in this case on the boundary of the interior ball H will 
be negative, owing to our choice of the normal z/*^. To construct an exam- 
ple with two boundary components, take hab = ^^^ab, where uj = 1 + 1/r; 
and let ^2 be the annulus Q = Br^ — B^^ where ro > 1 > ri > 0. Then 
we have that i? = 0, if = 2(1 + l/ro)~^(l - l/ro)/ro > on dBr^ and 
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H = 2{1 + l/ri)"^(l/ri — l)/ri > on 03^. Note that in this example the 
point i can not be at the origin r = 0. The generalization to k boundaries 
components is straightforward. 

Condition can be written in terms of physical quantities as K'^^PaVh > 
0. From equations (jSl)-® we see that, remarkably, every future marginally 
trapped surface on a maximal slice satisfies this condition. In contrast, con- 
dition (f29|l is a sufficient condition which is not expected to be also necessary 
in general. 

A future trapped surface is a time asymmetric concept. If we reverse 
the time direction —t°- we have — ^ —k"- and fc" — > —1°" and hence 

0+ — ©-, ©- Then, if we have a future trapped surface in one 

time direction, we will have 6+,B_ > in the opposite time direction and 
then it will be not trapped in this time direction. The singularity theorems 
[HK] implies that the development of a data which contain a future trapped 
surface will be geodesically incomplete into the future, but nothing is said 
about the past. In fact, in a realistic gravitational collapse, the spacetime is 
expected to have no singularities in the past. On the other hand, if the data 
have non-trivial topology (as the ones usually used in numerical simulations, 
see for example jlTj) we can apply the singularity theorem given in [22] 
to conclude that the spacetime will be geodesically incomplete in both time 
directions. This shows that the class of data constructed here is more general 
that the one with non-trivial topology, since it does not rule out data with 
regular past. 

It is perhaps surprising that 9_ and not 9+ is free data at the boundary. 
For example, we can not prescribe 9+ = and 9_ < at the boundary, since 
this will contradict inequality (jHUll . If we start with a boundary that satisfies 
9+, 9_ < 0, then, in general, a surface with 9+ = 0, 9_ < will be present 
in the interior of f2, this will correspond, for example, to the outer most 
trapped surface. But we can not use theorem 12.21 to start with such a surface 
at the boundary, the location of this surface can not be given a priori under 
our assumptions. On the other hand, theorem 12.21 allows to prescribe 9+ = 
at the boundary, this is just the time inversion of the case 9_ = 0, 9+ < 
which gives 9+ = 0, 9_ > 0. However, unless 9_ = 0, these surfaces will 
not be future marginally trapped. Moreover, these surfaces are located on 
the inner null boundary of the left quadrant of the Kruskal diagram (region 
IV of Fig. 6.9 in [iMj), and hence they are not expected to be present in any 
realistic gravitational collapse of matter. 

In the Kruskal diagram, trapped surfaces which satisfy condition pn|) are 
located in the left half of region II. Hence, they can be present in a spherically 
symmetric collapse. However, condition (n^Ull is a restriction in the sense that 
not in every collapse such surfaces will exist, as can be shown using the 
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spherically symmetric example. 

The choice of 9_ as a free data is dictated by the signs in the boundary 
condition (fT4| . The boundary condition should have the appropriate signs 
to use the maximum principle in order to prove that the conformal factor ip 
is positive (see sections El and Ej). 

Finally, we note that theorems 12.11 and 12.21 show how to obtain the whole 
solution in the exterior region f2 solving the equations in the compact domain 
f2. This can be used in numerical computations to solve for the whole data 
on a finite grid (see for example |2H1, [221) 

3 Preliminaries 

In this section we summarize some results from functional analysis and the 
theory of linear elliptic partial differential equations which we will use in the 
following proofs. Some of these results, although standard, are not easily 
available in the literature. 

The operators (defined after Eq. 0) and (defined by (|i8|) ) which 
appear in the Hamiltonian and Momentum constraints respectively, are lin- 
ear, second order, elliptic operators in divergence form. We will use the 
Hilbert space approach to the boundary value problem for this kind of op- 
erators. This approach has the advantage to be both simpler and applicable 
to a broader classes of solutions than other methods. 

Let fl a bounded open domain in M^, it will be assumed in the following 
that the boundary dQ is smooth. We shall use the following functions spaces 
(see [T], [21] for definitions, notations, and results) defined in fl: the set 
of m times continuously differentiable functions C"^{Q), the Holder space 
C""'"(fi), where < a < 1, the corresponding spaces C""(0), C™'°(0), the 
space C^{Q) of smooth function with compact support in Q, the Lebesgue 
space LP{Q), the Sobolev space if*(f2), and the local Sobolev space H^^^^Q) 
where s is a real number. 

An elliptic operator in divergence form has an associated bilinear form B. 

^ After this work was completed, there has appeared an article by D. Maxwell jHS] which 
studies the boundary condition 9+ = 0. In a new version of this article, which appeared 
after this work was submitted for publication, Maxwell makes an important improvement 
and he is now able to construct solutions which satisfy 8_ < 8+ = at the boundary. 
This is done imposing an extra condition which involves H and Kabv'^v''- Remarkably, the 
solutions obtained in ^| and the ones obtained in the present article do not, in general, 
overlap. However, it is important to emphasize that both set of solutions do not cover 
all possible black hole exterior regions. It is still an open problem how to construct and 
characterize all initial data for black holes exterior regions. 
I am grateful to D. Maxwell for useful discussions. 
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In the particular cases of and L/j these bilinear forms are given by (|37|l 
and (|l5|) respectively. A bilinear form is called weakly coercive over H^{Q) if 
there exist constants Ai > and Aq > such that following inequality holds 

^{y,v) > Ai||f ll^i(f^) - Ao||t;|||o(n) (31) 

for all V G H^iVL). If we can take Aq = in (p?T|l we say that B is strictly 
coercive. B is called bounded in H^iVt) if there exist a positive constant C 
such that 

\B{v,u)\ < (32) 
for all u,v E H^{Q). For a given B, we define the following spaces 

V = {ueX : B{v, u) = for all v G H\n)}, (33) 
W = {ueX : B{u, v) = for all v G H\n)}. (34) 

When B is symmetric we, of course, have V = W. 

The following theorem is the basic existence tool for the linear boundary 
problem, it is a consequence of the Lax-Milgram theorem and the Fredholm 
alternative in Hilbert spaces. 

Theorem 3.1. Let B be a bilinear form which is bounded and weakly coercive 
in if^(fi). Then V and W have both finite dimension. Moreover, let F 
be a bounded linear functional over H^{Q,). Then, there exist a function 
u G H^iVl) such that 

B{v, u) = F{v) for all v G H\n), (35) 

if and only if 

F{w) = 0, (36) 

for all w eW. 

We will apply this theorem for both second order elliptic equations and 
second order elliptic systems. Although the methods in both cases are, in 
many aspect, similar, second order elliptic equations have distinguished prop- 
erties, we will treat them separately. 

Consider the operator L^u. If we multiply by f G Cq{Q) and integrate 
by parts we obtain the symmetric bilinear form 

B(M,f)= / {h'^'^DauDhV - Ruv) dn, (37) 
Jq 

where dfi is the volume element with respect to hab- We have the following 
result, which will allow us to use theorem l3. II for the existence proof. 
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Theorem 3.2. Let B be given by ( T^TjJ . Assume that the metric satisfies 0. 
Then B is bounded and weakly coercive in H^{Q). 

For a proof see for example in [30| and [H] the same result is proved 
under more general assumptions on the coefficients. 

Consider the following linear boundary value problem for the elliptic op- 
erator Lh 

LhU = fmn, (38) 
NhU = (pon dn. (39) 

Note that our assumption on hab implies that H G C°°{dQ). A weak solution 
of this problem is a function u G H^{Q) which satisfies 

B{u,v)+ [ HuvdS = - [ fvdfx+ [ ipvdS, (40) 

for all V G H^{Q), where B is given by p7|) and dS is the surface element 
on dQ. One can check that a smooth weak solution u of (jiO|) will satisfies 
(f38| - (j39|) . This kind of boundary conditions are called natural boundary 
condition for the operator Lh, see [2], and fJW\ for a discussion on this 
method of treating general boundary problem. 
We define the bilinear form B' by 

B'{u,v) = B{u,v)+ HuvdS, (41) 
Jan 

and the linear functional F given by 

F{v) = - f fvdfx+ [ ipv dS. (42) 
Jn Jdn 

That is, B' and F are defined as the left and right hand side of Eq. (|40|l 
respectively. We want to apply the abstract existence theorem l3. II for B' and 
F. This theorem is not very useful unless we have a characterization of the 
null space W of B'. The maximum principle will be used to ensure that W 
is trivial under additional assumptions on the coefficients. 
We have the following version of the maximum principle. 

Theorem 3.3 (Weak Maximum Principle). Let B' given by ()1T|) . where 
the metric satisfies Assume that R > and H > 0. Let u G H^{Vt) 
satisfies B'(m, w) < (>)0 for all v G H^iSl), f > 0. Then u < (>)0 in Vt or 
u is a positive (negative) constant. 
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This theorem can be obtained from the more general results given in |37| . 
The proof is similar to the proof of theorem 8.1 in fSlj. 

We will also use the following standard versions of the maximum principle, 
see [211. 

Theorem 3.4 (Strong Maximum Principle). Assume that the metric 
satisfies ^ and that R > 0. Let u G H^iVt) satisfy LhU > in VL. Then, if 
for some ball B CC Q we have 

supw = sup-u > 0, 
B n 

the function u must be constant in Q. 

Theorem 3.5 (Hopf). Assume that the metric satisfies Q and that R> 0. 
Let u G C'^iyt). Suppose L^u > in Q. Let Xq G dQ be such that 

(i) u is continuous at Xq; 

(a) u{xq) > u{x) for all X E Q and u{xo) > 0; 
(Hi) d^l satisfies an interior sphere condition at Xq. 

Then the outer normal derivative of u at xq, if it exist, satisfies the strict 
inequality 

du , , 

- > 0, (43) 

We can now prove the basic existence result for the elliptic linear bound- 
ary value problem (|ifl|l . 

Theorem 3.6. Assume that the metric satisfies (jH). Assume also that R > 
0, H > and that either R or H is not identically zero. Let f G L'^{Q), 
(f G L'^{dVt). Then there exist a unique weak solution u G H^{Vt) of the 
boundary value problem ()4n|l . If, in addition, we have f < 0, ^ > 0, then 
u>0. 

Proof. Let B' and F be given by (|lT|l and (|i2|) . In order to apply theorem 
13. H we need to prove that F is bounded, B' coercive and bounded and finally 
that the corresponding null space W for B' is trivial. The boundedness of F 
follows directly from Holder inequality and the assumptions on / and (p. 

By theorem 13.21 we have that B is bounded, then to prove that also B' is 
bounded we only need to prove that the surface integral in ()1T|) is bounded 
in H^{Q). Using the generalized Holder inequality we obtain 



Huv dS 
an 



< 



\H\\L\dQ) I |m| |L4(an) I \v\\LHdn)- (44) 
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Since u,v E H^{Q), we can use the trace theorem (see (Tj) to conclude that 
u,v E //^/^((9f2), by the imbedding theorem we have that u,v E L^(9fi). 
Then, we can replace the L'^{dQ) norms in the left hand side of (|44|l by 
H^{Q) norms and hence the surface integral is bounded in H^{Q). 

Using that if > we have that B'(f,f) > B(t>,t>), then the coerciveness 
of B' follows directly from the coerciveness of B given by theorem 13.21 

From the maximum principle 13.31 it follows that every w E W should be 
constant, by the assumption that either i? or if is not identically zero it 
follows that w = 0. The last statement about positivity of u follows directly 
from the maximum principle 13.31 □ 

We note that if we make the stronger assumption H > we don't need 
to use the maximum principle in the previous proof, the uniqueness followed 
directly in this case. 

The following global estimates for the non homogeneous case will play a 
fundamental role in the non linear existence proof. 

Theorem 3.7. Assume that the metric satisfies (0). Let u E H^{Q) be a 
weak solution of (|^ . Assume also that f E L'^{fl) and ip E L°°{dfl). Then 
u E C"(n), < a< 1. 

The interior estimate part of this theorem (that is the fact that u E 
C"(n)) is standard. The crucial part is the regularity up to the boundary 
under the assumption ip E L°°{dfl). This can be proved using Campanato 
spaces, following similar arguments as in j23| and Theorem 3.8 of [25j, where 
the interior estimate is proved. 

We turn now to the elliptic system defined by Lhw". If we multiply by a 
E C^{fl) we get the following bilinear form defined over vectors 



We have the following analog to theorem 13.21 

Theorem 3.8. Assume that the metric satisfies 0. Then B defined by (|i5|) 
is bounded and weakly coercive on H^{Q). 

To prove this result, we first check that B satisfies the Legendre-Hadamard 
conditions and then we use the Garding inequality for functions with com- 
pact support. Using the assumption of the smoothness of dfl we can extend 
the result to H\n) using the Sobolev extension theorems. See for example 
[23] for a general discussion of this kind of bilinear forms. 

Finally, for the non linear existence proof we will use the following fixed 
point theorem (see for example [24]). 




(45) 
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Theorem 3.9 (Schauder fixed point). Let B be a closed convex set in a 
Banach space V and letT he a continuous mapping of B into itself such that 
the image T{B) is precompact, i.e. has compact closure in B. Then T has a 
fixed point. 

4 The Momentum Constraint 

The aim of this section is to prove theorem I2.1I The momentum constraint 
is solved using the York splitting as follows. Given a trace free tensor 
we look for a vector field w"" such that the tensor K"-'' defined by 

j^ab ^ gab _ ^C^^yb ^ (46) 

satisfies Then, it follows that the vector will satisfy the equation 

\.hWa = D'Qab in n, (47) 

where 

= D\Cw)ab. (48) 

We want to prescribe K°''^h'a^b freely at the boundary dfl, this suggests 
to impose the following boundary condition for equation (|48ll 

{Cw)aby'' = on dVL, (49) 

where Ch has been defined in (fT5ll . This boundary condition satisfies the 
Lopatinski-Shapiro conditions required in the elliptic estimates of [3]. More- 
over, this is a natural boundary condition for the operator L/i in the following 
sense. Assume that and are sufficiently smooth vector fields, we use 
(|i8|l and integrate by part to obtain the following identity 

/ UaLhw''d^l= [ iChwy'uaUbdS-B{u,w), (50) 
Jq J an 

where the bilinear form B(u, w) is given by From equation (fHUll we 

see that if L/^w" = in Q and {Cw)abJ^°' = on dil, then w°' is a conformal 
Killing vector. 

For the moment let us not consider the singular behavior at i. Then we 
have the following boundary value problem 

Lhw" = r in (51) 
{Chw)abv'' = on dVt. (52) 
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A vector field w G H^iVl) is a weak solution of the boundary value problem 
dnU-dni if 

B(m,w) = - / Ua.rd^i, (53) 
Jo, 

for all u" G H^{Vl). Consistently with (jSHI), we define a conformal Killine; 
vector field as weak solution of 

B(m,O = 0, (54) 

for all u°' G H^{Q). We have the following existence theorem. 

Theorem 4.1. Assume that the metric satisfies 0. Let J" G i^^(f2), p > 
6/5 he a vector field such that 

r^adfi = 0, (55) 



n 

for all conformal Killing vectors in f2. Then there exist a weak solution 

G H\n) of dp. 

Proof. We will use theorem 13.11 By theorem 13.81 we have that the bilinear 
form B defined by ()i5|l is bounded and coercive in H^{fl) under our assump- 
tions on the metric hat- Using Holder inequality and the Sobolev imbedding 
theorem one can prove that the functional 

F{u) = [ Uardfi, (56) 
Jn 

is bounded in H^{Q) if J" G ViVt), p > 6/5. Finally, condition (f55| is just 
equation F{^) = in theorem l3. II □ 

The boundary value problem (f5T|) - (f52|l has been studied in elasticity, usu- 
ally under stronger regularity assumptions on the coefficients (see for example 

ISH). 

We analyze now the singular behavior at the point i. As we point out 
before, it will be convenient to take advantage of the conformal invariance of 
the equations and to write them with respect to the metric (|18|1 . Since the 
Ricci tensor of the metric (fT8|l vanishes at i and we are in three dimensions, 
the Riemann tensor vanishes there too. Hence the connection and metric 
coefficients of (fTK| satisfy in normal coordinates x'^ (with respect to /i^^) 
centered at i 

r:^ = 0(r''), h',^=5,,+0{r"). (57) 
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In these coordinates let '^fiat ^ C°°{B^ \ {^}) be a trace free symmetric 
and divergence free tensor with respect to the flat metric 6ki, 

6.k^}l, = 0, d,^}l, = OmB,\{t}, (58) 

with 

^}L = 0(r'-^), dm, = 0{r'-'), dd^}l, = Oir'-'). (59) 

All these tensors have been characterized in theorem 14 of [Tn|. Denote by 
Qsinn the /I'-trace free tensor which is given by 



Qsl, = Xe i^flat - I K, ^flat h'^')- (60) 



The tensor Qgi^g has two important properties. The first one is 



D'.Q±, = 0{r'-'). (61) 

This is a consequence of equations (|H7|l and The second is the following. 
If the metric hat admit a conformal Killing vector ^"'j then 

/ ^bDaQ:lg dfi = P''ka + rS, + Aa + (P^ "(0 + Q-) g„ (62) 
Jn 

where the constants k"-, S"", a and are the conformal Killing data at i for 
given by (fT6|) and the constants P", J°, A and Q"^ partly characterize the 
tensor field Qg^^g (there exist two free functions that can be freely prescribed 
in Qsing, see [IB] for details). Equation (f62|l is written in terms of the metric 
hab, it is possible to write it in term to the rescaled metric /i^^, in this case 

For the regular part Q^lg we will assume that 

Qf,geC^{Q\{z}), (63) 

and 

gab ^ oir-^), dQ"'' = Oir~^), ddQ"^ = 0{r~^). (64) 

An example is of course Q'^lg G C°°{fl), however the most general condition 
matches naturally with the fall off of Qsing in the sense that the constant 

gives the fall of 0(r~^), the constants A and J"^ the 0(r~^) one and the 
constants P" the 0{r~*), see [T6] . 

Using the fall off (|6i|l we obtain 

ibDaQtgdli= I QtgibyadS-\im f Q^gibUadS, (65) 
Jan JdB^ 

Qtg^^yadS. (66) 

an 
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We can now prove theorem 12. II 
Proof of theorem 12.11 

Proof. Since the conformal factor Uq defined by (fTTjl is smooth, we can write 
the equation with respect to the metric /i^^ 

J,^,w"' = r in fi, (67) 
{Ch,w')aby'^ = on dVL. (68) 

Define J" = D'^^Q^K By the ^ and ^ we have that J" G ^^(fi), with 
6/5<p<3/2. Then we can apply theorem l4.11 In the presence of conformal 
Killing vectors the condition F{^) = in theorem l4. II can be written as the 
equation (pT| using equations (jH^ and fHHll . This shows the existence of 

The regularity (i.e. w;" G C°^{n\{i}) and consequently K'''' e C°^{Cl\{i}) 
) follows directly from the standard elliptic regularity theorems of 
Remains to show that K"-'^ satisfies 

IK"''] = 0(r-^). (69) 

We have that 

li^'^'l < ig^'l + K^'w;)"")!. (70) 

The first term on the right hand side of this inequality satisfies, by (fHn| 
and (fHij) . the estimate (|69|l . We only need to estimate the second term. Let 
fl' GG Q, we have the following inequality 

sup \{Ch'wy'')\ < C\\w\\ci(n>\B,) < C\\w\\h3{Q'\b,), (71) 
xen'\B^ 

where the first inequality is obvious and the second is a consequence of the 
Sobolev imbedding theorem. Using the interior elliptic regularity theorem of 
in] we have 

\\w\\H3(n'\B,) < C{\\J\\H^n'\B,) + \ \'w\\l2{q'\b,))- (72) 

Since w G H^{Q) then | |w| |2,2(q'\bo is clearly bounded for all e. For the first 
term in the right hand side of (f72|l we use the assumptions f59l) and f64ll to 
conclude that 

C 

\\J\\m{n'\B,) < ^, (73) 
for every e. Combining the previous inequalities it follows that 

sup \{CH'wr')\<C/r\ (74) 

xGn'\Bt 

for every e. Then, using (f7n|l the desired result follows. □ 
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5 The Hamiltonian Constraint 



In this section we will prove theorem 12 .21 We begin with the item i) of this 
theorem, that is the existence of the function ipo defined as a solution of the 
boundary value problem (f22|) - (f2^ . This function will play an important role 
in the non linear existence proof; it is essentially the Green function of the 
operator with the boundary condition (f23|) . 

In the following theorem 6i will denote the Dirac delta distribution with 
source at i. 

Theorem 5.1. Assume that the conformal metric hat satisfies ^ and that 
dQ is smooth. Assume also that R > in Q, H > on dfl and that either 
R or H is not identically zero. Then, there exist a unique solution ipo of the 
boundary value problem (f22|) - pij) . Moreover, ipo satisfies 

Lhipo = -ATiSi, (75) 

in Q, ipo > in Cl \ {i}, ipo G C°°{Cl \ {i}) and it has the following form 
V'o = Xt/^ + u with u G H'^iyi). 

Proof. Observing that 1/r defines a fundamental solution to the flat Lapla- 
cian, we obtain 

A(^) = -47r5, + x, (76) 

where x is a smooth function on VL with support in \ -Be/2. The ansatz 
^0 = Xth' + ^ translates the original equations into the following equations 
for u 

LhU = -Lh{^) - X in fi, (77) 
NhU = on dVL, (78) 

where we have deflned Lh by the expansion Lh = ^ + Lh in normal coor- 
dinates centered at i. A direct calculation shows that Lhixe"^'^) £ L'^i^t) fl 

{pL\{i}) . Then, by theorem l3.6[ there exists a unique solution u G H^{VL) 
of (f77| - (f78|) . We can use the standard elliptic regularity theorems to con- 
clude that u G H'^{fl), which in particular implies (by the Sobolev imbedding 
theorem) that u G C^{fl). We can further use the elliptic regularity in the 
region fl \ {i} to conclude that u G C°°{Q \ {i}). 

To show that ipo is strictly positive, we observe that it is positive near i 
(because r~^ is positive and u is bounded). Take e small enough such that 
ipo is positive on dB^^. Note that Lhtl^o = and ipo is smooth in fl\B^. Let 
Xq E Q\Bf^ the point where ipo reach its minimum in Cl\B^. Since Lhipo = 0, 
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we can use the strong maximum principle 13.41 to conclude that either ipo is 
constant in ri\B^ or Xq is in the boundary ofri\B^ and il^lx) > ipi^Xo) for all 
X E Q\B^. Since we have assumed that either R or H is not identically zero, 
ipo can not be constant. Assume that ipo{xo) < 0. Since > on dB^, then 
it follows that xq G dQ. Consider the function —ipo- We will apply theorem 
13.51 on ri\B^ for this function to get a contradiction. Hypothesis i), ii) and 
iii) of this theorem are satisfied. Using the boundary condition (f78|l we have 

Av''Dai-i')ixo) = -H{-ij)ixo) < 0, (79) 

which contradict (j^ . □ 

We note that theorem l3.5l is essential to prove the strict inequality ipQ > 0, 
if we use theorem 13.31 we only get ipo >0. 

We treat now the general, non linear, case. Set = ipo + u, then the 
boundary value problem given by (fTI|l and (fT^ can be written as follows 

LhU = f{x,u) in Q, (80) 
NhU = ip{x,u) on dQ, (81) 

where 
and 

cp(x,u) = -{iJo + uf^i + r-j — r~v5' (S3) 

where we have used the notation 

e_ = -ipi, KabV'^v' = (^2, (84) 

to emphasize that the functions ipi and in the following are arbitrary, 
non-negative functions which are not necessarily given by by (IMj) . 

We first recall some special properties of the functions /(x, u) and ip{x, u). 
The function f{x,u) can be written as 

/(x, Ml) - f{x, U2) = (mi - U2)f{x, Ml, M2), (85) 

where ^ 

f{x,u,,U2) = ^^^J2^iJo + u,y'\i^o + U2)-'-'. (86) 

j=0 

Clearly f{x,Ui,U2) > for any Mi,M2 > 0. Also, for any m > we have 

0>/(a;,n)>/(x,0). (87) 
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For (p{x, u) we have a similar expression 

ip{x, Ml) - ip{x, U2) = - Ui)(f{x, Ml, M2); (88) 

where 

2 

<^(x,Mi,M2) = ^V5i(V'0 + Uif'^^lpo + + ^P2ilp0 + Ml)"^"^(V^O + M2)~^"^- 

i=o 

(89) 

If we assume ipi, ip2 > 0, then we have that 0{x, Mi, M2) > for any Mi, M2 > 0. 
However, in contrast to f{x,0), the function ip{x,0) has, in principle, no 
definite sign. 

We define m+ as a solution of the following linear problem 

LhU+ = fix, 0) in Q, (90) 
NhU+ = ip2i^o^ on on. (91) 

li ip2 > and f{x,0) G L'^{fl), then, by theorem I3.6[ there exist a unique 
solution M+ > of (|nni)-(inH). Note that the function ^1 defined by (l25ll-(IT7ll 
is given by ipi = ipo + u^. The following theorem will prove items ii) and iii) 
of theorem 12.21 

Theorem 5.2. Assume that the hypothesis of theorem 15. 1\ are in force. 
In addition, we assume that y^i,y^2 o,re smooth function on d^l which sat- 
isfy (pi,(p2 > 0, KabK"''^ipQ'^ G L^(r2) and (p{x,u+) > where m+ is the 
unique solution of ((001) ^(jHIJ ■ Then, there exist a unique positive solution 
u e H'^{Vl) nC^lnX {i}) of the boundary value problem M-CT. More- 
over, < M < M+. 

Proof. The proof is based on the Schauder fixed point theorem 13.91 Define 

B C CO(fi) as 

B = {ue C°(fi) : < M < M+}. (92) 

Clearly B is closed and convex subset of the Banach space C^ipL). Let w E B, 
we define the map T{w) = m as follows. Let u be the unique solution given 
by theorem 13. 6| of the following linear boundary value problem 

LhU = f{x, w) in Q (93) 
NhU = ip{x,w) on dn. (94) 

Using equation ()88|1 . ()85|1 . the maximum principle 13.31 and the condition 
<^(x, M+) > one shows that T{B) C B. Using the elliptic estimate 13.71 we 
have that u G C^i^Q), and hence the image of T{B) is precompact because 
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C"{fl) is compactly imbedded in 0^(0,). It is also clear that T is continuous. 
Then, by the Schauder fixed point theorem 13 .91 there exist a fixed point 
T{u) = u with u G C^ipL). The standard elliptic regularity implies that in 
fact u G H^{Q) n C^ifl \ {i}). This finishes the existence part. 

To prove uniqueness, let assume that there exist two solutions ui and U2, 
the difference Ui — U2 will satisfy the equations 

Lh{ui - U2) = {ui - U2)fix, ui, U2) (95) 

Nh{ui - U2) = {U2 - Ui)ip{x,Ui,U2) (96) 

where / and (p are given by (j^ and fHnil respectively. Assume that ui 7^ U2 
in some set. Since ui and U2 are continuous functions, we can find a set 
Q' G fl such that Ui > U2 in VL' (we change Ui by U2 if necessary) , Ui — M2 = 
on r, where V denotes the portion of dVL' contained in dVt and U1—U2 > 
on r. By (j95ll we have Lh{ui—U2) > 0. U1—U2 can not be a positive constant 
in Vt' because / > and > 0. Then by the strong maximum principle there 
exist xq G dO! such that {ui — U2){xq) > {ui — U2){x) for all x G fl'. The 
point Xq can not be on dfl' — T because there we have Ui — U2 = and this 
will contradict Mi > U2. Then Xq G F. Using the boundary condition (|96|l we 
get 

Au'^Daiui - U2)ixo) < -H{ui - U2){X^) < 0. (97) 
We use theorem 13.51 to get a contradiction. □ 

Note that ip{x,Uj^) > is precisely condition (^ni) . To finish the proof 
of the item iii) of theorem 12.21 it remains only to prove the inequality (f3n|) . 
Using equations (fT2|l - (fT3|) we obtain 

0_ - e+ = 2Nh^ = 2if{x, u) > 2if{x, u+) > 0. (98) 

Using the sub and super solution method, boundary value problem of the 
form has been studied in under stronger regularity assumptions 

on the coefficients (see also [34]). In [2I| a related equation has been also 
studied using a variational approach. 
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